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We discuss theoretically quantum interface between light and a spin polarized ensemble of atoms 
with the spin > 1 based on an off-resonant Raman scattering. We present the spectral theory of 
the light-atoms interaction and show how particular spectral modes of quantum light couple to 
spatial modes of the extended atomic ensemble. We show how this interaction can be used for 
quantum memory storage and retrieval and for deterministic entanglement protocols. The proposed 
protocols are attractive due to their simplicity since they involve just a single pass of light through 
atoms without the need for elaborate pulse shaping or quantum feedback. As a practically relevant 
example we consider the interaction of a light pulse with hyperfine components of Di line of ^'^Rb. 
The quality of the proposed protocols is verified via analytical and numerical analysis. 
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I. INTRODUCTION 

Light-matter quantum interface is a basic element of 
any quantum information network aiming at long dis- 
tance quantum communication, cryptography protocols 
or quantum computing, see ^ (2| ■ Light is a natural car- 
rier of quantum information and a macroscopic atomic 
system can be efficiently used for its storage. Critical in- 
gredients of the quantum interface are the quantum mem- 
ory and entanglement protocols, which allow high-fidelity 
interchange (transfer, storage and readout) of quantum 
states between the light and relatively long-lived atomic 
subsystems. A number of promising theoretical propos- 
als for the high-fidelity memory and entanglement pro- 
tocols has been put forward, which can be classified as 
based on off-resonant interaction, such as Raman inter- 
action 3, quantum non-demolition (QND) interaction 
with quantum feedback - QNDF [4], multiple QND in- 
teractions Ref. Q , and on resonant interaction using elec- 
tromagnetically induced transparency - EIT [H, [7|- In 
spite of a large number of proposals, experimental real- 
ization of the complete storage plus retrieval quantum 
memory with the fidelity higher than classical is yet to 
be achieved. High fidelity storage (but not complete re- 
trieval) has been demonstrated via the QNDF approach 
and as light-to-atoms teleportation Q. Low fidelity 
storage and retrieval of the coherent and single photon 
pulses based on EIT 10] and Raman [ll| processes have 
been recently demonstrated. 

All experiments on light-atoms quantum interface to 
date are conducted with alkali atoms. Previously the 
off-resonant quantum interface protocols Q utilized 
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light with the detuning higher than the hyperfine split- 
ting of the excited state. Under this condition magnetic 
sublevels of a hyperfine state of an alkali atom can be ef- 
fectively reduced to a spin-half system. As shown in the 
present paper eliminating this restriction and thus us- 
ing the complete magnetic multipole system of the alkali 
atom ground state opens up new possibilities for quan- 
tum memory and entanglement generation with ensem- 
bles of such atoms. The quantum description of corre- 
lations in coherent Raman scatte ring in the Heisenberg 
formalism has a long history [l^ . [isl . 0, UBl • However 
in quantum information applications the Raman process 
is often discussed in terms of formal A-scheme config- 
uration. In contrast in our treatment we describe this 
process in terms of polarization multipoles for multilevel 
ground state alkali atoms: gyrotropy (orientation), lin- 
ear birefringence (alignment) and higher multipole com- 
ponents. The field subsystem is described by a set of 
polarization Stokes components. This approach allows 
for complete description of the mean values and fluctu- 
ations of the light and optically thick atomic medium. 
Although examples treated in this paper concern atoms 
initially pumped to one level corresponding to a single A- 
scheme (as shown in figures [T^nd (2]) , our formalism also 
allows to treat more complex initial states consisting of 
several coupled A-schemes. The formalism is also ap- 
plicable beyond the Heisenberg equations in a linearized 
form. Finally, the multipole expansion facilitates system- 
atic and compact calculation of the coupling constants in 
the effective Hamiltonian taking into account the full hy- 
perfine structure of the relevant atomic levels [T6| . 

We shall discuss the Raman quantum memory scheme 
first considered in 0]. We shall add the retrieval step to 
the protocol and analyze the complete procedure for a 
realistic model, taking full account of the multilevel hy- 
perfine and magnetic sublevel structure of an alkali atom. 
We derive the polarization-sensitive coupled wave-type 
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dynamics of an optically thick atomic ensemble and op- 
tical field. Using spectral mode decomposition for light 
and spatial mode picture for the atomic ensemble we 
show that the quantum states of light and the quasi-spin 
of atoms can be effectively swapped or entangled. The 
general mathematical formalism for such atomic system 
with angular momentum higher or equal than one, has 
been developed in Refs.[l^[l3|. An important advantage 
of the memory and entanglement schemes in such a sce- 
nario is in that they can be realized in a single pass of 
light through atoms and without any feedback channel, 
which makes experimental realizations more feasible. 

As an elementary carrier of the quantum information 
we consider a squeezed state of light. The relevance of 
quantum memory for these states for quantum informa- 
tion processing has been illustrated in proposals [l8] . We 
introduce Stokes operators for a light pulse consisting of a 
circularly polarized strong classical mode and the quan- 
tum squeezed light in the orthogonal circularly polar- 
ized mode. We show that for the quantum Stokes vari- 
ables there is a convenient symmetric interaction with 
the alignment tensor components of atoms, whose spin 
angular momenta are equal or greater than one. In the 
quantum memory scenario the quantum information of 
the light subsystem can be effectively mapped into the 
alignment subsystem of atoms. In the entanglement sce- 
nario the excitation of atomic spins with coherent pulse 
generates a parametric-type interaction process, which 
results in creation of the strong correlations between the 
quantum fluctuations of the Stokes components of the 
transmitted light and the alignment components in the 
spin subsystem of atoms. 

We consider the interaction via the Di line of ^^Rb, as 
an important example, where a convenient spin-one sys- 
tem exists in the lower hyperfine sublevel of its ground 
state. For verification of the memory protocol we discuss 
its figure of merit for the mapping of the input squeezed 
state and for its retrieval after reading the quantum copy 
out from the spin subsystem with a second coherent light 
pulse. We discuss how the fidelity for the proposed quan- 
tum memory channel could be defined and compared 
with the respective classical benchmark based on direct 
measurement of the squeezed state parameters. We show 
that in an optimal configuration the quantum fidelity is 
always higher than the limit for the optimal competing 
classical channel. 

Throughout our analysis we neglect coupling of the 
ground state atomic degrees of freedom to other variables 
than the forward propagating light field. This approxi- 
mation is justified on time scales which ar short com- 
pared to dephasing time of the ground state coherence. 
Neglecting coupling to light modes propagating in other 
than the forward direction can only be a good approxi- 
mation in the limit of high optical depth as detailed in 
section UlIAl Finally, the influence of the atomic motion 
capable of washing out spatial spin patterns is not taken 
into account, which restricts our analysis to cold atomic 
samples. 



II. OVERVIEW OF THE PHYSICAL 
PROCESSES AND BASIC ASSUMPTIONS 

A. The schemes of experiments under discussion 

In this paper we will consider two alternative experi- 
mental situations respectively shown in figures [T] and [2l 
For both schemes a 100% right-hand circular polarized 
classical light pulse interacts with an ensemble of ultra- 
cold spin-oriented atoms. The principal difference is in 
the direction of the collective spin of the atomic ensem- 
ble, which can be oriented either along (Fig[T]) or opposite 
(Figl2) to the propagation direction of light. 

Consider first the process shown in figure [TJ For an 
off-resonant right-hand polarized pulse, such that inco- 
herent scattering losses do not frustrate the spin polar- 
ization, there will be no interaction between the light and 
atoms. However if a small portion of a left-hand polar- 
ized quantum informative light prepared in an unknown 
squeezed vacuum state is admixed to the classical coher- 
ent pulse, the coherent scattering channel will be open. 
The portion of the weak quantum light will be coherently 
scattered into the strong classical mode. In this coherent 
Raman process the polarization quantum subsystem of 
the probe light and the spin subsystem of atoms can ef- 
fectively swap their quantum states. The quantum state 
can be mapped into the alignment-type fluctuations of 
the spin subsystem and further stored in the form of a 
certain standing spin wave for relatively long time. It can 
be readout on demand with a second probe light pulse. 

In the experimental situation shown in figure [2l for 
a transparent medium a small portion of the left-hand 
polarized photons will emerge as a result of the coher- 
ent elastic Raman scattering of the strong light pulse. 
This process creates entanglement the left-hand polarized 
quantum modes of light and the alignment- type quantum 
fiuctuations in the atomic spin subsystem. This kind of 
deterministic entanglement can be interesting for imple- 
menting a quantum repeater protocol between remote 
atomic systems. 

The processes shown in figures [T] and [5] have practi- 
cally identical theoretical description. Thus we will dis- 
cuss them in parallel paying attention to the differences 
in their physical nature. We begin by introducing the 
relevant set of the field and atomic variables involved in 
these processes. 



B. Field subsystem 

The polarization states of the light subsystem can be 
specified in terms of the following flux-type variables for 
different Stokes components. The Heisenberg operators 
for the total photon flux considered at the spatial point 
z and at the time moment t are given by 

%{z,t) = -^E(-)(z,t)E(+)(z,i) (2.1) 
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FIG. 1; Schematic diagram showing the geometry of the pro- 
posed experiment on quantum memory and the scheme of 
relevant excitation transitions, see text for details. 



where E^^)(z, t) are the positive/negative frequency com- 
ponents of the electric field Heiseberg operators. We as- 
sume the quasi-monochromatic and forward propagating 
probe light such that the spectral envelope of the modes 
has a carrier frequency (D, and the light beam has a cross 
area Sq. 

The Stokes components responsible for three alterna- 
tive polarization types are given by 



Si(z,t) 



S2(z,t) 



S3(z,t) 



Sqc 
27:% u) 



Sqc 



E^^\z,t)E^^\z,t) 
-Ei-\z,t)Ei+\z,t) 
'_Ei-Hz,t)Ei+Hz,t) 

(2.2) 
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These components subsequently define imbalance in the 
photon fluxes with respect to the x/y Cartesian basis 
(S3), to the ^/r) Cartesian basis rotated at 7r/4-angle 
(Si ) and to the R/L basis of circular polarizations (S2), 
see [g. 

The Stokes variables (|2.2p obey the following commu- 
tation relations 



FIG. 2; Schematic diagram showing the geometry of the pro- 
posed experiment on entanglement between light and atomic 
subsystems and the scheme of relevant excitation transitions, 
see text for details. 



where eijk = ±1 depending on the order of the indices 
i 7^ j 7^ k. The (5-function on the right-hand side of this 
commutation relation is an approximation valid if only a 
finite bandwidth of the field continuum is important for 
the correct description of the low frequency fluctuation 
spectrum. This assumption is in accord with the rotating 
wave approximation which we will further use in the de- 
scription of the atom-field interaction. For the processes 
shown in figures [T] and [2] the commutation relations (|2.3[) 
can be simplified to 



S3(z,t),Si(^',t) = 2icS{z 



(2.4) 



Si(z,t),^j(z',t) 



2ieijkc6{z - z')^k{z,t) (2.3) 



where for the case of small fluctuations the operator on 
the right-hand side is replaced by its expectation value. 
This Stokes component is approximately conserved in the 
linear evolution, and because of the conservation law for 
the number of photons in the coherent process one has 
S2 = Sq. The component Sq is the integral of motion 
in our model and Np = SqT, where T is the interaction 
time, gives the average number of photons participating 
in the process. 



C. Atomic spin subsystem 

The angular momentum polarization of a-th atom can 
be conveniently described in the formalism of the irre- 
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ducible tensor operators 



2K 

2F, 



^ E «kqI^oM')(Fc 



M',M 



KC 



2Fo + l 



FoMKQ ^KQ 

(2.5) 



With these expansions the set of atomic dyadic operators 
originally defined in the subspace of the Zeeman states 
|FoM), where Fq is the total (spin + nuclear) angular 
momentum of atom and M is its Zeeman projection, is 
transformed into the set of the relevant irreducible tensor 
operators. HereC " are the Clebsh-Gordan coefficients, 
with the rank and the projection Q which respectively 
vary in the intervals: < K < 2Fo and —K < Q < K, 
see Ref.0. 

The Heisenberg dynamics of the irreducible compo- 
nents (j2.5p preserves the following commutation relations 



T, 



X V[l 



K" 



I'o I'o ro 



X (-1)2^0+^'' c^S::Q,f^")^,,(t) (2.6) 

For the interaction processes shown in figures [T] and [2] 
only the following two combinations of the alignment op- 
erators and the longitudinal orientation component con- 
tribute to the dynamics of the atomic spin fluctuations 
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where Fz°'\t) is the Heisenberg operator of the angular 
momentum projection for a-th atom on z-axis. For this 
set of operators the commutation relation (|2.6p leads 



= SablCiF^'''>{t)+SablC3f^o\t) 



5ab 2if^;\t) 



(2.8) 



The coefficients ci and C3 are given by 
3 



Cl 



C3 



Fo(Fo + l)(2i^o + l) 

Q[{Fo~\){Fo 



2)] 



1/2 



[7^^o(^^o + l)(2Fo - l)(2i^o + l)(2i^o + 3)]'/' 

(2.9) 



and the higher-rank irreducible component T30 weighted 
with coefficient C3 contributes only for atoms with Fo > 
1. 

The operators Txy^ and T^"^^ have clear physical mean- 
ing and define the components of alignment tensor with 
respect to either x, y or ^,77 Cartesian frames, which 
were earlier introduced with definition of the polariza- 
tion Stokes components of light. These operators can 
also be related with the transverse components of a quasi- 
spin, which could be defined in the system of two Zeeman 
states coupled via the A-type excitation channels shown 
in figures [T] and [21 However in context of our discussion 
it is more important to emphasize the tensor nature of 
these components, which are responsible for the fluctu- 
ations of linear birefringence of the sample initiated by 
either stimulated or spontaneous Raman processes shown 
in these flgures. 

The important step in further description of the col- 
lective dynamics of the atomic spins consists in assump- 
tion that the original discrete distribution of the point- 
like atomic spins can be smoothed with procedure of the 
mesoscopic averaging. Then for a mesoscopically thin 
layer located between z and z -I- Az one can deflne 



Z<Za<Z + Az 



Z<Za<Z + Az 



where the sum is extended on the atoms located inside 
the layer. These expressions deflne the smoothed spatial 
distribution of the alignment components for which the 
first line of commutation relation (12.811 transforms to 



%y{z,t),f^rj{z',t) ^ S{z - z')iCi3Tz 
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Fo(Fo-f l)(2Fo + l)(2Fo-f 3) 



Cl3 



(2.11) 



where in the case of small fluctuation the right-hand side 
is replaced by an expectation value for the density of 
atomic angular momentum J-"^. As in the previous sec- 
tion we assume that this projection of the collective spin 
is approximately conserved in the Raman process such 
that the total angular momentum F^, = — ±FoNa, 
where L is the sample length and N^^ is the number of 
atoms. Since the incoherent losses are ignored Na re- 
mains the integral of motion in this model. 



D. Dynamic equations 

The master equation governing the dynamics of atoms- 
fleld variables in the Heisenberg picture can be derived 
similarly to the way it was done in Ref. [lB|- We apply 
the effective Hamiltonian derived in that paper and ne- 
glect the dissipation channels caused by incoherent scat- 
tering. 
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We will assume that the squeezed probe radiation for 
the excitation scheme, shown in figure [U is characterized 
by a degeneracy parameter (i.e. by a mean number of 
photons in the coherence volume of the quantum radia- 
tion) higher than unity. As was shown in Rcf. 21] there 
is a special scenario when the A-system is probed by ex- 
tremely weak SPDC light source consisting of strongly 
correlated photon pairs, which should be cooperatively 
scattered by all atoms of the ensemble. This process 
cannot be described by the Hamiltonian of Ref. [l^| and 
should be discussed separately. 

Omitting the derivation details we present the follow- 
ing Heisenberg wave-type equations describing the tem- 
poral and spatial dynamics of the Stokes components and 
the alignment components 



Ei{z,t) = KiS3(z,t) - 2eE2%y{z,t) 
'E.3{z,t) = -KiSi(z,t) -t- 2e'~2%^{z,t) 

%y{z,t) = ~Clf^ri{z,t) + ci3.eTz'E.i{z,t) 

T^jj{z,t) ^ VL%y{z,t) - ci3eTz'B3{z,t) 

(2.12) 



" d 


1 d 


dz 


cdt 


' d 
dz 


1 d 


cdt 



d_ 

di 
d_ 
di 



These equations are written for the physical obsrevables 
and have a rather clear physical structure. They show 
that in a linearized regime the modification of the light 
and atomic polarization in a forward passage results in 
a combine action of gyrotropy and linear birefringence. 
In the discussed case the former manifests itself as an 
averaged interaction between atoms and light and the 
latter actually creates the quantum fiuctuation interface 
between these subsystems. Let us point out that in a 
formal description based on a A-type configuration the 
fluctuations of birefringence (alignment) components are 
usually treated as an atomic coherence created by Raman 
process in a quasi-spin subsystem and the importance 
of the phase-matching gyrotropy can be underestimated. 
We show how this effect can be properly taken into con- 
sideration by direct solution of Eqs. (|2.12p . 

Equations (|2.12p introduce the following important pa- 
rameters: Ki, n = 2Qq + Oi and e. The first two de- 
scribe the changes in the mean classical parameters of 
light and atoms, ki describes the gyrotropy of the sam- 
ple caused by the average angular momentum orientation 
of the atomic spins, fli is the light shift between the Zee- 
man sublevels Fq, M = ±Fo and Fq, M = ±Fo =F 2. It is 
additionally assumed that the atoms interact with an ex- 
ternal magnetic field and JIq is the relevant Zeeman split- 
ting for the neighboring sublevels. The third parameter e 
is the most important and is responsible for the coupled 
dynamics of the fiuctuations of the polarization compo- 
nents of light and of the atomic alignment components. 
The microscopic expressions for all these characteristics 
for the case of hyperfine transitions of alkali atoms are 
summarized in Appendix El 



Before we discuss the solutions of these equations we 
turn again to the physics of the processes shown in fig- 
ures [T] and [51 The difference between these processes is 
only indicated by different signs of J^z in the second pair 
of equations (|2.12[) . For positive (FiglT]) these equa- 
tions yield the normal wave-type solution leading to the 
swapping of the quantum states between the light and 
spin subsystems. Then the memory protocol can develop 
in the following way. If the quantum field is prepared in 
an unknown left-hand polarized squeezed-vacuum state 
its momentum/position quadrature components can be 
expressed as 



Xp cx -ie'^a + ie-^'^a^ 



(2.13) 



where the phase 9 determining the orientation of the 
squeezed and anti-squeezed quadratures in the phase 
plane as well as their variances are the unknown parame- 
ters of the state. If such an unknown squeezed state is su- 
perimposed with a classical coherent mode one straight- 
forwardly gets that 



li oc Xp 

'-3 OC Xq 



(2.14) 



The crucial point for the atoms-field dynamics described 
by equations (|2.12p and for the memory protocol is that 
the direction of the axes x, y and ^, 77 remains completely 
unknown in the experiment, see figure [T] Thus the un- 
known squeezed state transforms into the unknown polar- 
ization squeezing such that the Si component is squeezed 
and S3 is anti-squeezed. The dynamics of these compo- 
nents is further developing in accordance with equations 
(|2.12p and under certain conditions the squeezed state 
of light can be converted into the squeezed state of the 
alignment components of the atomic spin angular mo- 
menta. 

For the negative sign of J-'z (Pig^ there is no nor- 
mal wave solution of the equations (|2.12p . As we will 
show in this case there is an exponential enhancement 
of the quantum fluctuations of the Stokes components 
and of the atomic alignment. Such an anomalous spin 
polariton wave describes an entangled state of these sub- 
systems. After interaction the low frequency modes of 
the polarization and intensity spectrum of the outgoing 
light will be entangled with the alignment-type standing 
wave modes of the atomic spins. Such process can be for 
example utilized in the quantum repeater protocol with 
remote atomic systems. 



III. 



COUPLED DYNAMICS OF THE ATOMIC 
AND FIELD SUBSYSTEMS 



The Heisenberg equations (|2.12p can be solved via the 
methodjjf Laplace transform similar to how it was done 
in Ref. 16|. The solution and veriflcation of the commu- 
tation relations are given in Appendix |B] In this section 



6 



we discuss the physical consequences which are mainly 
important from the quantum information point of view. 
For the sake of simplicity and without loss of generality, 
see Ref.[l6'], we will ignore the retardation effects and 
restrict our discussion to the low frequency fluctuations 
such that in equations ()2.12p we neglect time derivatives 
of the field variables. We also restrict our consideration 
to the special case when the relevant pair of Zeeman sub- 
levels, which is either |Fo, Fq) and |Fo, Fq — 2) (Fig[T|) or 
|Fo, — i^o) and \Fo,—Fq + 2) (FigH]), is always degenerate 
such that Cl = 2r2o + ili = during the interaction and 
Oo = without the interaction. This non-critical but 
convenient simplification allows us to consider storage of 
quantum states of light in the time independent stand- 
ing spin wave without any Zeeman-type oscillations, see 
appendix |B] 

The qualitative character of the cooperative dynamics 
of the atomic and field variables can be illustrated by the 
dispersion relation (|B6[) for the spin polariton modes and 
light modes defined in the semi-infinite medium and for 
infinite interaction time. Then temporal and spatial dy- 
namics of the process can be expressed by the respective 
Laplace modes s and p introduced by expansion (|B4p . 
The Laplace variables for the inverse transform can be 
parameterized as s = —iQ,p = iq. Then the inverse 
transform reveals a spectral Fourier expansion over the 
set of relevant temporal and spatial modes. The wave- 
type dynamics makes these modes coupled via the fol- 
lowing dispersion law 



A 
q 



A = ~2ci3e'S2J^^ 



(3.1) 



This dispersion relation reflects the coupling of low fre- 
quency temporal fluctuations of light to fine scale spatial 
fiuctuations of atoms and of fast temporal fluctuations of 
light to long scale fluctuations of atoms respectively. For 
an atomic medium of the length L and a light pulse of 
the duration T, the combination ATL, characterizing the 
overall coupling strength between light pulse and atomic 
ensemble, quantifies this scale relation in natural units. 



A. Quantum memory protocol 

In the following we apply the input/output relations 
derived in appendix [B] to identify the suitable parame- 
ters for writing and retrieving to and from a quantum 
memory. For similar systems the optimization of a quan- 
tum memory has been analyzed very thoroughly in the 
single and few- mode situations f?, Ts'l , while here we are 
interested in a multi-mode situation considering pulses of 
squeezed light of a duration much longer than the inverse 
bandwidth of squeezing. Prior studies of squeezed state 
storage in [l^ considered a steady state regime, while 
here we discuss a temporal sequence of writing, storage 
and retrieval. 



1. Write-in stage 

Consider the experimental scheme shown in figure[T]for 
a pulse-type excitation of the system with a portion of 
quantum light with the following correlation properties 



1 _lz 

e "1 ) ^2 



1 
2^ 



(3.2) 



where {...,.. .}_|_ denotes the anti-commutator of two ob- 
servables. These correlation functions describe the out- 
put generated by the intra-cavity subthreshold degener- 
ate parametric amplifier, see Ref.^^l- The freely prop- 
agating light with such correlation properties can be as- 
sociated with a squeezed state of a harmonic oscillator 
where Si component is squeezed (^ < 0) and S3 compo- 
nent is anti-squeezed ($3 > 0) such that (1 + Ci)(l + C3) = 
1 . The spectral bandwidth of the informative part of the 
quantum light is limited by the longest of two correlation 
times 



Tl 
T3 



IC 

2 

ic 
2 



(3.3) 



where 7c is the cavity loss rate through the output mir- 
ror and is the efficiency of the downconversion pro- 
cess. For high level of squeezing one has T3 ^ ti. This 
inequality indicates that the minimal duration T of the 
quantum light pulse should be considerably longer than 
the longest time T3 = Tc. 

The general solution (|B8p . given in the appendix [B] 
can be straightforwardly rewritten in the following form 

Ti°^^{z) = cosKiz7^j,(z,T) + sin Kiz '7^,,(z, T) = 
ciseT^J^ dtJo[2[-A{T-t)z]'/^) Si"(t) . . . 

7J^"'(z) = cos Ki 27^^(2, T) — sin kiz 'I'xy{z, T) = 
-ci^eT, j^dt Jo(2[-A(r - t)z]i/2) ^^(t) + . . . 

(3.4) 

which shows how the input field operators are transferred 
into space-dependent atomic spin operators. The dots in 
the right hand side stand for the contributions of the 
atomic operators responsible for reproduction of input 
atomic state, which make this transformation imperfect. 

The principle condition, which makes the memory pro- 
tocol feasible, is that for the low-frequency spatial fluc- 
tuations on the left hand side of Eas. (|3.4p the contri- 
bution of the input spin fluctuations on the right hand 
side is suppressed. In the spectral expansion of the op- 
erators 7j°"*(z) and 7Jf"*(z) such fluctuations belong to 
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the spectral domain q qc = [^-—[^]^^^- In accordance 
with the dispersion law p.ip . which is asymptotically 
valid for T,L oo, the transform (|3.4|) performs map- 
ping of the field fluctuations within the spectral interval 



The noise contribution caused by the vacuum fluctua- 
tion of the input fleld can be suppressed in the low fre- 
quency domain ft < fl'^ ^ [l^y^Ji/s^ see Eq.dBTj. At 
the retrieval stage of the protocol the spin fluctuations 
< n < onto the spatially dependent spin fluctua- within the spectral interval < g < are mapped onto 



tions within the spatial spectral domain oo > g > |^|tc. 
For reliable mapping of the field quantum state onto the 
atomic alignment subsystem the latter domain should es- 
sentially overlap with the spectral domain where the in- 
put spin fiuctuations are suppressed. For the efficient 
memory protocol the following inequalities should be ful- 
filled 



\A\tc < qc 



L ^ -t: qc 



\A\T 
L 



1/2 



(3.5) 



These inequalities imply that the contribution of the in- 
put spin fluctuations are suppressed and ensure that the 
spatial fluctuations of 7^°"*(z) and 7^"*(z) correctly re- 
produce the temporal dynamics of the input fleld fluctu- 
ations (j3.2[) in the relevant part of their spectrum. How- 
ever the dispersion relation (|3.ip has a singular behavior 
for the function 17 = il{q) near the point q such that 
the shot noise part of the input fluctuation spectrum is 
reproduced at the zero point of spatial spectrum. Thus 
near the spectral point q — > 0, associated with an integral 
collective mode of the standing spin waves in the limit 
L — > oo, the informative part of the fluctuation spec- 
trum will be invisible. However, for the sample with a 
finite length L the squeezed state can be mapped onto 
the integral collective modes of 7J°"'(z) and 7^"*(z) if 



\A\tc < L^' < qc 
quantum light. 



i.e. for a spectrally broad incoming 



2. Retrieval stage 

For retrieval of the quantum state back onto light the 
atomic ensemble should be probed with another strong 
coherent light pulse. Optimal retrieval occurs when 



KiL = 2n X any integer 



(3.6) 



the time dependent field fiuctuations in the frequency 
spectral domain oo > Q ^ l^'l^c, where the correlation 
length is given by Ic ^ q^^ ^ L. At the retrieval stage 
of the memory protocol the following inequalities should 
be fulfilled 



\A'\h « n[ 



\A'\L 



T' 



1/2 



(3.8) 



which, after the obvious replacement of temporal and 
spatial parameters, are symmetric to the inequalities 
p.Sp and have similar physical meaning. Thus the time 
dynamics of the output field operators S°^* (t) and S™* (t) 
correctly reproduce the spatial distribution of spin fiuc- 
tuations in the relevant part of the spectrum. 



3. Physical requirements and numerical simulations 

Basic requirements for the proposed memory proto- 
col can be summarized by combining the inequalities 
p.Sp and (|3.8p where parameters A and A' are given by 
Eq. (|3.1[) . The following limitation should be imposed on 
the number of atoms Na and on the numbers of photons 
in the coherent strong field Np and N'p participating in 
the process at the write-in and retrieval stages respec- 
tively 



c'^NaNp > 1 



(3.9) 

The substitution of Ic from (|3.5p into p.Sp yields 

N'p < Np (3.10) 

From evaluation of incoherent losses we obtain the fol- 
lowing conditions on the numbers of atoms and photons 



Then following (|B7p the output field operators are given 

by 

Sr (0 - -2e^'^j\z Jo(2[-A'(L - z)tY'^) fr\z) + ... ^^ere a 



7Vpf7+ _2, N'pa+^_^ « So 



(3.11) 



(3.7) 



where the initial state of the atoms is '7^™*(z) and 7Jg"*(z) 
modified in the write-in stage of the protocol. The dots 
on the right hand side again indicate the imperfection of 
the transform because of the presence of the input field 
operators associated with the vacuum modes. 



is the cross section of the incoherent scattering 
for a left-hand polarized light from the atoms populating 
the Zeeman sublevel M = Fq and crpg_2 ^^"^ cross sec- 
tion of the incoherent scattering for a right-hand polar- 
ized light from the atoms populating the M = Fo — 2 Zee- 
man sublevel; 5*0 is the beam/sample cross area. The first 
line indicates that the atomic medium should be trans- 
parent for the informative quantum radiation. The sec- 
ond line indicates that the incoherent scattering should 
have negligible infiuence on the dynamics of the spin co- 
herence during both the write-in and the retrieval cycles. 
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These requirements lead to the following basic demand 
on the density of atoms np, resonance radiation wave- 
length A and the sample length L: noX^L :$> 1, i.e. the 
on resonance optical depth of the sample needs to be 
large 

Optimization of the inequalities (|3.10p and (|3.1ip re- 
quires that Np < Na < Np. To explain this choice 
let us again consider the coherent Raman scattering as 
the physical background for the memory protocol. In the 
write-in stage of the protocol the most desirable process is 
the annihilation of the photons of the quantum modes as 
a result of the coherent scattering into the classical mode 
followed by the transfer of their quantum state onto the 
alignment of atoms. The swapping process can be easier 
done for the photons in the wings of the spectral distribu- 
tion (|3.2p and is more difficult to realize for the resonance 
photons. That is why the process runs more efficiently 
with a stronger coherent field with Np > Na- At the 
retrieval stage the coherent Raman scattering takes pho- 
tons from the strong mode into the quantum modes. The 
spin fluctuations of the atomic alignment tensor transfer 
back into the outgoing modes of the quantum light. To 
make this process more efficient, i.e. to increase the num- 
ber of scattering events, it is desirable that the sample is 
extended and contains a large number of the atoms such 
that Na> N'p. 

In figure [3] we show how the spectral variances of the 
polarization components of light and atoms are modi- 
fied after sending the broadband squeezed light (tc 0) 
through the atomic sample. The graphs clearly display 
the swapping mechanism at the write-in stage of the 
memory protocol. The graphs are normalized to the vac- 
uum state variance and the deviations from it arc ex- 
pressed by the respective Mandel parameters 1 + for 
i = 1, 3, which depend on either frequency (for light) 
or wave number q (for atoms). The input squeezed state 
is described by two spectrally independent parameters 
1 = 10 and 1 ^i'^ = 0.1. The suppression of the 

correlations in the low frequency domain of the temporal 
spectrum of light is compensated by the enhancement of 
the correlations in the atoms for the low frequency part 
of the spatial spectrum. The dominant role of the collec- 
tive modes in this process is a result of the approximation 

Tc^O. 

The data, presented in these graphs, corresponds to 
the red detuning of light near the Z?i-line of ^^Rb for 
atoms in the hyperfine level _Fo = 1 a-nd for two selected 
values of the cooperative interaction parameter ATL — 
— 10,-40. The first value corresponds to the detuning 
—205 MHz from Fq — 1 ^ F = 1 resonance, where ki = 
(no gyrotropy effect) and for approximately ten percent 
losses estimated from Eqs. p.lip . The second number is 
achievable for the detuning of a few thousand MHz in the 
red wing of the Fq ~ 1 ^ F ~ 1 transition and for the 
same level of losses. 

Figure [3] shows the spectral variances of the polar- 
ization components for the atoms and light at the re- 
trieval stage of the protocol. The figure shows how well 




<n 



0.00 3.14 6.28 9.42 12.57 15.71 

nT 




in 



0.00 3.14 6.28 9.42 12.57 15.71 

FIG. 3: The spectral variances of the Stokes components (up- 
per panel) and of the atomic alignment components (lower 
panel) before and after the interaction with the broad- 
band squeezed light for cooperative parameter ATL = — 10 
(squares) and ATL = —40 (triangles). The black solid lines 
indicate the original fluctuation spectra in the light and spin 
subsystems. The gray lines in each of the panels indicate the 
original spectra for a complementary system. 



the recovered state can reproduce the input. The pa- 
rameters are chosen such that A'T'L = —2 correspond- 
ing to ATL = -10 and A'T'L = -8 corresponding to 
ATL — —40. Let us point out once more that for the best 
overall efficiency, the number of photons in the strong co- 
herent pulse at the retrieval stage should be smaller than 
the number of photons applied at the write-in stage of 
the protocol and smaller than the number of atoms. As 
follows from these results the retrieved quantum state of 
light can reproduce the input state only in certain parts 
of the fluctuation spectrum. Further optimization should 
allow for identification of the best temporal mode, where 
the retrieval of the original squeezed state would be op- 
timal. 

Figures [H] and [5] show how the spectra of the light and 
atoms are modified at the write-in and retrieval stages 
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FIG. 4: The spectral variances of the alignment components 
of the atoms (upper panel) and of the retrieved Stokes com- 
ponents of the light (lower panel) after the readout of the 
state stored in the atoms at the write-in stage of the mem- 
ory protocol, see Fig|3] The values of cooperative parameter 
A'T'L = —2 (squares) and A'T'L = -8 (triangles) are coor- 
dinated with the data of FigO The black solid lines indicate 
the original fluctuation spectra in the spin and light subsys- 
tems. The gray lines in each of the panels indicate the original 
spectra for the complementary system. 



FIG. 5: Same as in Fig|3] but for the finite-bandwidth 
squeezed light. The ratio of interaction time to the corre- 
lation time of squeezing for these graphs is T/tc = 10. 



difficult to store and retrieve than the broadband input 
state, the fact already established for the storage phase in 
3] . This is the direct consequence of the imperfection of 
the swapping mechanism for mapping the low frequency 
fluctuations, as follows from the dispersion law discussed 
above. 



with the finite-bandwidth squeezed light as the input. 
The dependencies plotted in these graphs illustrate the 
importance of the dispersion relation between the tempo- 
ral and spatial modes participating in the protocol. The 
ratio of the interaction time to the correlation time of 
the squeezing is taken to be T/tc = 10. As follows from 
the displayed results, for the finite-bandwidth squeezed 
light the integral collective modes are not optimal in both 
storage and retrieval steps of the protocol. A certain op- 
timization procedure is necessary for the best encoding 
of the quantum information into particular domains and 
modes of the fluctuation spectra for the light and atoms. 
Another important observation which follows from these 
graphs is that the finite bandwidth squeezing is more 



4-. Fidelity 

Fidelity is the fundamental criterion for any quantum 
memory or teleportation scenario. It can be relatively 
simply defined for a single mode situation and for a pure 
input state, and is more subtle for a general case. In the 
following we show qualitatively that in the present case 
the quantum scheme has always a better fidelity than 
for what can be called a classical memory /measurement 
protocol. 

We assume that the classical memory scheme utilizes a 
balance homodyne detection of the input squeezed state 



10 



IS 

> 



<u 
a. 

0) 




0.00 



3.14 



6.28 



9.42 



12.57 



15.71 



qL 




15.71 



FIG. 6: Same as in Fig|4] but for retrieval of the finite- 
bandwidtli squeezed state mapped onto tlie spin subsystem, 
see FigO The original spectrum for the squeezed light in the 
lower panel relates to the case of T' = T. 



for the measurement of its parameters. The measurement 
should yield the degree of squeezing and the directions of 
the linear polarizations shown in figure [TJ Such a proce- 
dure is displayed in figure [7^ as an overlap of the Wigner 
function associated with the input state with its mea- 
sured counterpart. This overlap in the limit of high level 
of squeezing can be expressed by the following fidelity 



F: 



1 



v/1 + (D^eN)^ 



(3.12) 



where D3 — (1 + ^3)/2 is half of the variance for the 
anti-squeezed S3 component, which is assumed to be re- 
liably defined, and On tt/N is the angular uncertainty 
remaining after N measuring attempts to identify direc- 
tion of x,y or ^,77 axes. Following Ref.jl^l the following 
constraint should be fulfilled in order to make the classi- 
cal protocol efficient 



< [DsONf < 1 



(3.13) 




(b) 



FIG. 7: The overlap of the Wigner function for the input 
squeezed state (green ellipse) with its "classical" (a) and 
"quantum" (b) counterparts (red ellipses). 



where T/v ^ T/N is the duration of each measuring at- 
tempt as a fraction of the total measuring time T. It is 
clear that for D3 —^ 00 the number of measurements N 
should grow up to infinity such that for any limited time 
T the constraint (|3.13p can never be fulfilled. 

On the contrary, in the quantum memory protocol it 
is only necessary to effectively reproduce the variances 
associated with the squeezed and anti-squeezed polariza- 
tion components. The directions of the x, y and ^, 77 axes 
remain completely unknown in the interaction process. 
This is shown in figure [7)3 with a more effective overlap 
than in case of figure [7^. The fidelity can be conveniently 
written in terms of Mandel parameter for the input and 
output states 



F : 



[(2 + ei" + er)(2 + ei"+e3™*)] 



1/2 



(3.14) 



where "out" characteristics should be associated with 
optimally defined spatial mode. Comparing l|3.12p and 
p.l4p for the same ratio T/tc and for the high level of 
squeezing one can expect that in the optimized case the 
fidelity for the quantum protocol will be always higher 
than the relevant classical benchmark. 



B. Atoms-field entanglement 

The entanglement between the light and the atomic 
alignment subsystems created in the process shown in fig- 
ure [2] is attained in a result of the transformations (|B12p . 
(IB13P under the constraint (|3.6p . The entanglement can 
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be written in the following canonic form 



Jq 
^ h{t)ErHt)dt+ I g{z)f^^\z)dz 





(3.15) 



The alignment spin waves "/^""'(z) and 7^"*(z) are de- 
fined by the upper lines of Eqs. (|3.4p . The standard entan- 
glement between the quadrature components of certain 
collective field and atomic spin modes can be recognized 
in Eqs. (|XT5)) , see also relations (P?T^ and 

Formally the temporal mode h{t) and the spatial mode 
g{z) can be found by solving the following integral equa- 
tions 



h{t) 



dt' 



' AL ' 


1/2 






t' -t 





h [2[AL{t' -tf"^] h{t') 



ciaeJ^. / dzh{2[A{T-t)zf'^] g[z) ^ Q 



9{z) 



r dz' 


AT 




' z 


z' ~ z 



h(2 [AT{z'-zf 
o(2[A(L-z)t]i/') h(t)^ 



9iz') 
(3.16) 



Due to the presence of non-symmetric integral operators 
in the upper lines of these equations their actual solution 
does not responsibly exist in terms of real functions. One 
can search for a pair of real functions h(t) and g{z) which 
minimize the absolute value of the left hand side of (|3.16p . 
This circumstance is reflected by the " limit to zero" in the 
right hand side, which can be approximately approached 
for certain optimal observation conditions. 

It is intuitively clear that this can be achieved for a sys- 
tem consisting of a large number of atoms {e\Tz\L ^ 1) 
and of a large number of coherent photons {eE.2T ^ 1). 
Using the common properties of the Fredholm-type equa- 
tions the solutions for the spatial and temporal modes of 
p.l5p can be found in terms of the eigenfunctions of the 
integral equations p.l6p . 



solution under rather general assumptions. This allows 
for a convenient and clear interpretation of the Heisen- 
berg dynamics of atoms and field at every step of the 
discussed quantum information protocols. 

We consider two basic protocols: quantum memory 
and quantum entanglement. The memory protocol is an 
example of the swapping mechanism between the atomic 
and field subsystems. The quantum information, which is 
originally encoded in the polarization degrees of freedom 
of the light wave, can be mapped onto the spin stand- 
ing wave associated with atomic alignment components. 
However because of the multimode nature of the inter- 
action process the quantum information can be spread 
among a number of spatial spectral modes. We show 
how the protocol, particularly at the final retrieval stage, 
can be optimized. As an important step of such optimiza- 
tion we have comprehensively discussed the imperfection 
of the quantum memory channel and made qualitative 
comparison of its figure of merit with a competing clas- 
sical memory /measurement scheme. 

The second protocol is generation of entanglement in 
the spin oriented atomic ensemble probed with a cir- 
cularly polarized coherent light mode. Forward scat- 
tered Raman photons appear strongly correlated with an 
alignment-type coherence of atomic scatterers. When the 
process becomes extended in space as well as in time it 
creates pairs of entangled temporal and spatial modes. 
We introduced the system of two integral equations, 
whose solution defines the structure of these modes. 
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IV. CONCLUSION 



We discuss the polarization sensitive interaction be- 
tween the Stokes components of light and the alignment 
components of the atomic ensemble as a resource for 
quantum information interface. Our model accurately 
describes the quantum nature of interaction between the 
spin subsystem of ultracold alkali atoms and a plain light 
wave propagating through the sample, in the absence 
of losses. The method has two important advantages. 
First, the process does not require any special manipula- 
tions like a feedback or adjustment of the driving classi- 
cal light. This creates good outlooks for its experimental 
implementation. Second, the model yields an analytical 



APPENDIX A: COUPLING PARAMETERS OF 
THE HEISENBERG EQUATIONS 

The wave-type Heisenberg equations (|2.12p straight- 
forwardly follow from the commutation relations between 
the operators of the respective quantum observables with 
an effective Hamiltonian responsible for the process of 
coherent forward scattering. The latter was derived in 
Ref . . The exact equations are linearized with as- 
sumption of small fluctuations when S2 and Tz are ap- 
proximately conserved and interpreted as classical values. 
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Then the gyrotropy constant is given by 

F - - I 



[Fo(^o + l)(2Fo + l)]'/' 



(Al) 



Here o'p'^piOj) is the dimensionless orientation income 
into the polarizabihty tensor for the Fq ^ F hyperfine 
transition 



ap^^p{u) _ (-1) ^\f, F, 



111 

F 



SqC -h{u! - LOfFo) 



(A2) 



where dpoF is the reduced dipole moment and lofFo is 
the transition frequency. The hght shift is given by 



n -(1) f-\ V3^2 

ih = > ap p(uj) 77:7 

V ° '[Fo{Fo + l){2Fo + l)f' 



(A3) 



This expression is similar to definition of the gyrotropy 
constant (|Aip because both the parameters come from 
the same Faraday-type interaction term of the effective 
Hamihonian. 

The couphng constant e responsible for the alignment- 
type interaction is given by 



(A4) 



(2) 

where ap^p{uj) is the dimensionless alignment contribu- 
tion into the polarizabihty tensor 



\d 



1 1 2 

Fq Fq f 

2 



FoF\ 



Attuj 



SqC -h{uj - LUpPg) 



(A5) 



All the parameters are the spectrally dependent values 
such that Ki = Ki((D), = Qi{lu), e — e((D) and this 
dependence can be very important in practical calcula- 
tions. 



APPENDIX B: THE SOLUTION OF THE 
HEISENBERG EQUATIONS 

We derive the solution for those physical conditions 
when only low frequency temporal fluctuations are con- 
sidered as a quantum information carrier. That lets us 
completely ignore the retardation effects associated with 
a finite sample size. Practically this means that we can 
neglect the time derivations in the first two lines of the 
system ((2?T2t . 

As a first step we make the following local rotational 
transformations for the Stokes components of the field 
subsystem 

Si(z,i) = cos (^(z, t) Si(z, <) — sin(p(z, i) ^3(2, i) 
S]ii(z,i) = sin(p{z,t)Ei{z,t) + coa(p{z,t)E3{z,t) (Bl) 



and for the alignment components of the atomic subsys- 
tem 

Ti{z,t) = cos(p(z, t)%y{z,t) + sini^(z, t)T^j^{z, t) 
Tts_{z, t) = - sin Lp{z, t)%y{z, t) + cos (p{z, t)T^ri{z, t) (B2) 

where (f{z, t) = K,iz+Vtt. Without losing of generality the 
parameter Cl can be set as to ignore any non-principle 
freely precession of atomic spins associated with external 
magnetic field and light shift. It is also convenient to set 
the sample length L as kiL — {2tt x any integer), then 
the pair of Stokes variables Si(z, t), 'B.^^z^ t) and Si(z, i), 
53(2, t) will coincide at the output of the sample. 

Applying the orthogonal transformations (|B1[) and 
(IB2[) to the system (|2.12l) the latter can be rewritten in 
the following form 

^Si(z,i) = -2e%fi{z,t) 
oz 

d . - - 

— 5n(2,i) = 2e^2%K{z.,t) 
oz 



d_ 

dt 
d 



Ti{z,t) = ci3eTz^i{z,t) 



dt 



Tw{z,t) = -ci3eJ^z^m{z,t) (B3) 



Then its solution can be found by the method of Laplace 
transform. 

To show this we can define the Laplace images of the 
space-time dependent Stokes components of the probe 
light and of the collective alignment components of atoms 



5i(p,s) 



00 />oo 



JO 

oc poo 



^0 



dz dte-P^-"^ %{z,t) 
dzdte-P''-''^%{z,t) (B4) 



with i = I, H and /i = I, H. Then equations (jBSp is 
transformed to the system of algebraic equations, which 
in turn can be straightforwardly solved similar to how 
it was done in Appendix C of Ref.[ll]. The important 
parameter for this procedure is the determinant of the 
system, which is given by 



A(p,s) = [sp + 2ci3e2S2J^.]' 



and its pole 



A 
P 



(B5) 



(B6) 



with A = —2ci3e^'E,2^z, can be linked with the dispersion 
law for the spin polariton wave propagating through the 
sample, see Eg. p.ip and discussion around it. 

Finally for A < one obtains the following solution for 
the Stokes components of the light subsystem 
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EiiL,t) 





'-AL 


fdt' 




/o 


t-t' 



11/2 



.h[2[-AL{t-t')f'^) 



2e^2 / dz J( 



',(2[-A{L~z)t]"^)j^%z) 





'-AL' 


1/2 


fdt' 




JO 


t~t' 





Ji(2[-AL(t-0]'^')2S(0 

l/2> 



+2eS2 y dzJo[2 [-A{L - z)tY'' ) T^\z) 



(B7) 



and for the alignment components of the spin subsystem 



7i(z,r) = friz) - / dz' 



' -AT' 


1/2 






z — z' 





.'m1/2 



+ci3e^. / Jo ( 2 hA(r - t)zf^ ) 



tni(z,r) = fi'iz)- / dz' 



" -AT' 


1/2 






z — z' 





Ji(2hAT(z-z')]'^') 

1/2^ 



-CiseTz dtJo[2 [-AiT - t)zY'' ) S'^{t) 



(B8) 



Here L is the sample length and T is the interaction 
time. The solution performs an integral transform of the 
Heisenberg operators for the Stokes components in the in- 
cident light and of initial Schrodinger operators for 
the alignment components of atomic spins T^iz). The 
kernels of the transform are expressed by the cylindrical 
Bessel functions of the zeroth Jo(. . .) and first Ji(. . .) or- 
der. The original form of atomic and field operators can 
be recovered with the orthogonal transformations reverse 
to jHH) and (fB2ll . 

The derived solution preserves the following commuta- 
tion relations 



SM(z,i),Si(z,i')J = 2i5{t-t')% 
%^{z,t),fi{z',t)\ = -icxz5{z- z')T, (B9) 



The commutation relation for the Stokes components, 
given by the first line, differs from the original commu- 
tation relation in form (|2.4p . That is direct consequence 
of our ignoring the retardation effects. The argument of 



(5-function, which generally is t — t' — {z — z')/c, can be 
only approximately reproduced in the assumptions we 
did. The solution also obeys the following important 
cross-type commutation relations between the field and 
atomic variables 

Si(z,t),'7^(z',i') = -2ici3egip52-^z 

X Jo (2 [-A{z - z'){t - t')]'^') m - t') - e{z' - z)] 

(BIO) 

where the matrix g^^^ (with i, /i = I, H) is given by 



0, 1 

1, 



(Bll) 



These commutation relations clear indicate that the 
atomic and field variables always commute before and 
after interaction. 

For A > one obtains the following solution for the 
Stokes components of the light subsystem 
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1/2 


fdt' 




Jo 


t~t' 





/i 2[AL(t-t')] 



1/2 







1/2 






'0 







,1/2 



and for the alignment components of the spin subsystem 



f,{z,T) ^ friz) 



dz' 



AT 


1/2 






z ~ z' 





h[2[AT{z-z')f^)fr{z') 



+ci3e^, j dtio (2 [A{T - i)^]'^') 
1/2 

%{z,T) = ftiz 



dz' 



h(2[AT{z~z')f')U-{z') 



-ciae^. dtIo(2 [A(T - t)zf^) 



(B12) 



(B13) 



Here the functions /o(- • ■) and ■ •) arc the modified 
Bessel functions of the zeroth and first order. The solu- 
tion obeys the commutation relations (|B9p and the cross- 
type commutation relations are now given by 



^i{z,t),Tfj,{z',t') = ~2ici3e gif,^2^z 
x/o (2 [A{z - z'){t - t')]^^^) W - t') - 0{z' - z)] 



(B14) 
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